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ABSTRACT 

PSR J1740— 3052 is a young pulsar in orbit around a companion that is 
most likely a B-type main-sequence star. Since its discovery more than a 
decade ago, data have been taken at several frequencies with instruments 
at the Green Bank, Parkes, Lovell, and Westerbork telescopes. We measure 
scattering timescales in the pulse profiles and dispersion measure changes 
as a function of binary orbital phase and present evidence that both of 
these vary as would be expected due to a wind from the companion star. 
Using pulse arrival times that have been corrected for the observed peri- 
odic dispersion measure changes, we find a timing solution spanning 1997 
November to 2011 March. This includes measurements of the advance of 
periastron, w, and the change in the projected semimajor axis of the or- 
bit, i, and sets constraints on the orbital geometry. From these constraints, 
we estimate that the pulsar received a kick of at least ^50 km/s at birth. 
A quasi-periodic signal is present in the timing residuals with a period of 
2.2 X the binary orbital period. The origin of this signal is unclear. 

Key words: binaries: general - stars: early-type - stars: mass-loss - pul- 
sars: general - pulsars: individual: PSR J1740— 3052. 



1 INTRODUCTION 

Since their discovery more than 40 years ago, pulsars 
have proven to be incredibly useful and precise astro- 
physical laboratories, due in large part to our ability to 
time their exceptionally steady rotations. In particular. 



E-mail: madsense@phas.ubc.ca 



by timing a pulsar in orbit with a binary companion, we 
learn not only about the pulsar, but about the companion 
and how it affects the pulsar's orbit and radiation. This 
provides a unique method for probing properties of both 
the pulsar and the companion. 

Of the '--^ITO radio pulsars known to be members of 
binary systems, so far only four have been found whose 
companion is a main-sequence star. PSR B 1259— 63 or- 
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bits a type B2e companion |johnstonet"ai]|l992l). and 
PSR J0045-7319 has a Bl V companion dKaspi et alj 
Il994h. The discov ery of PSR J1638-4725 was reported 
bv lLorimer et alJ I2OO6) and was a result of the Parlce s 
multibeam survey for pulsars ( [Manchester et alJIlOOlh : 
its companion spectral type has no t yet been determin ed. 
PSR J1740-3052, discovered bv lStairs etall j200lh in 
the same survey, orbits a companion that has been iden- 
tified as most likely a late O-type or early B-type star 
Istairs et al.ll200ll : lBassa et alj|20ll1) . A coincident late- 
type star has be en shown not to be the companion 
( lTametalj|2O10h . The pulsar has an orbital period of 
231 days, spin period 570 ms, and characteristic age 0.35 
Myr. 

In this paper, we present an updated timing model 
for PSR J1740— 3052 using more than a decade of data 
collected at several observatories. These observations 
are outlined in Section (2] along with the updated tim- 
ing model and an analysis of scattering and dispersion 
caused by a wind from the companion. We see that 
dispersion-measure (DM) variations agree very well with 
those predicted by a simple stellar wind velocity model. 
In Section[3]we discuss secular changes in orbital param- 
eters and the presence of a low-frequency, quasi-periodic 
signal in the timing residuals, and we conclude in Section 
111 



2 DATA REDUCTION AND TIMING 
2.1 Observations 

Data were collected between 1997 November and 2011 
March on the Lovell Telescope at Jodrell Bank, UK, 
the Parkes telescope in New South Wales, Australia, the 
NRAO Green Bank Telescope in West Virginia, USA, 
and the Westerbork Synthesis Radio Telescope (WSRT) 
in the Netherlands, as listed in Table[T]along with the cor- 
responding total bandwidths, channel bandwidths, and 
observing frequencies. Data from Parkes and Jodrell 
Bank up to 2000 Nove mber 23 were used in the timing 
solution of lStairs et al] (2001). 

At Jodrell Bank, observations are made using ei- 
ther an analogue filterbank (AFB) or a digital filterbank 
(DFB). The AFB systems summed the orthogonal polari- 
sation channels in hardware and used I-bit digitisers. The 
DFB system uses 8-bit digitisers and generates folded 
pulse profiles with 1024 bins per pulse period and full 
Stokes-parameter information in real time. Total inten- 
sity profiles were formed in off-line analysis. 

Observations at Parkes are also made with either 



Table 1. Observing frequencies and bandwidths of the 
J1740— 3052 data. Where multiple values appear in more than 
one tield on a single line, they may be matched in the order 
listed. 



Frequency 


Total bw 


Channel bw Typical obs. 


(MHz) 


(MHz) 


( MHz ) time ( min ) 


Jodrell Bank AFB (1997 Nov - 2010 Apr) 


1376 


64, 96 


3 


24 


1380 


64, 96 


3 


60 


1396 


64 


1 


24 


1402 


64 


1 


24 


Jodrell Bank DFB (2009 Jan - 2011 Mar) 


1373.875 


128 


0.25 


24 


1381.5 


112.75 


0.25 


24 


1524 


384, 512 


0.5 


24 




Parkes AFB (1998 Oct - 2009 Mar) 




660 


32 


0.125 


25 


1374 


288 


3 


10 


1390 


256 


0.5 


10 


1518 


576 


3 


10 


2350 


288 


3 


25 




Parkes DFB (2005 Dec - 2010 Aug) 




1369 


256 


0.25, 0.5 


10 


1433 


256 


0.5 


10 


Westerbork Synthesis Radio Telescope (1999 Oct- 


2010 Feb) 


1375 


80 


0.156 


30 


1380 


80 


0.156 


20 


Green Bank BCPM (2001 Sept - 2009 Jan ) 


570 


48 


0.5 


20 


575 


48 


0.5 


20 


590 


48 


0.5 


20 


820 


48 


0.5 


20 


1190 


48 


0.5 


10 


1400 


96 


1 


10 


1660 


48 


0.5 


15 


1780 


48 


0.5 


15 


2200 


96 


1 


15 



an AFB or one of several DFB systems. The sampling 
time on the Parkes AFB is 250 |j.s, and the DFB pro- 
files had 1024 bins per period. The PDFB data were flux- 
calibrated using a pulsed noise diode signal injected into 

the feed. 

The BCPM ("Backe r et aljl997l) at Green Bank was 
a filterbank with flexible channel widths and 4-bit sam- 
pling. Profiles observed using the BCPM have a sam- 
pling time of 72 \is and are flux-calibrated using a pulsed 
noise diode. 

WSRT data are sampled every 409.6 |i.s using the 
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PuMa backend dVoute et ai] l2002h and are not flux- 
calibrated. 

Data obtained at every site are dedispersed as- 
suming constant DM at values between 739 and 
743 cm"'' pc. These variations in assumed DM do not 
affect the quality of the data products used here. 

Individual pulses are folded at the predicted 
topocentric pulse period to produce pulse profiles with 
high signal-to-noise. Folding is performed online, prior 
to calibration, in the DFB cases, and offline, after cali- 
bration, for the AFBs, BCPM, and WSRT. 

Except in the case of the DFBs, data from two 
orthogonal polarisations are summed — during the fold- 
ing process for the WSRT and the BCPM, and directly 
through the hardware for the AFBs. The DFB polarisa- 
tions are summed to form total-intensity profiles before 
obtaining pulse times-of-arrival (TO As). 



2.3 Timing analysis 

TOAs are processed using the pulsar timing soft- 
ware TEMPoQ with barycentric corrections applied us- 
ing the Jet Propulsion Laboratory DE405 solar system 
ephemeris |standish 1998) and clock corrections applied 
using the UTC(NIST) timescale. The uncertainties asso- 
ciated with the TOAs input into TEMPO tend to be under- 
estimates, and thus we increase the error on sets of data 
from different observatories and in different frequency 
ranges such that each set has a reduced of ~1. The la 
errors on parameters output by TEMPO are also doubled, 
providing error estimates we consider to be conservative. 

FITSCATTER-produced TOAs are used for the 
lower-frequency data which have scattering tails that are 
quite long, as discussed in Section [Z2l but we find that it 
makes no significant difference to the parameter values 
or reduced of the TEMPO fit (prior to error-scaling) 
for TOAs with frequency >1 GHz and only serves to 
reduce their correlation with higher-frequency residuals 
that do not show scattering. For this reason, we do not 
use FITSCATTER TOAs for data at frequencies > 1 GHz, 
which are the TOAs used in the full timing solution. 



2.2 Scattering and arrival-time fitting 

To obtain TOAs, a Gaussian standard profile is first 
constructed from a fit to high signal-to-noise data from 
each different instrument, and for the various frequency 
ranges on each ins trument, using the program BEIT 
( lKrameretal]|l994l) . Data with minimal scattering are 
used where possible, and BEIT is used to fit out the effects 
of scattering if necessary. Each pulse profile is cross- 
correlated with the appropriate standard profile to pro- 
duce a topocentric TOA. 

Data observed at low frequencies show significant 
scattering, particularly at frequencies <1 GHz. This is 
seen as a decaying exponential tail convolved with th e 
pulse profile, as seen in Figure 1 of IStairs et alj ( 1200 ih . 
The program EITSCATTER (written by MK) can be used 
to read in scattered pulse profiles along with a Gaussian 
standard profile and output both TOAs that account for 
the shifted peak induced by the convolution and scatter- 
ing timescales. 

We believe the scattering timescales found by the 
fits to be reasonable estimates. However, the TOAs pro- 
duced using EITSCATTER are less reliable than TOAs that 
do not account for scattering, and are not ultimately used 
in the full timing solution, as discussed in the next sec- 
tion. 



2.4 Orbital dependence of DM and scattering 
timescale 

Variations in disp ersion measure (D M) near periastron 
were detected by IStairs et al] 1 I2OOI1) and taken as evi- 
dence for a wind from a nondegenerate comp anion. The 
wind velocity model of iKaspi et al.l { 1996bh and refer- 
ences therein is used there to support the idea that the 
companion is an early B-type star rather than a late-type 
supergiant. Variations in DM in this model can be ex- 
pressed as 

where M is the mass-loss rate of the companion in 
Mq yr~^ and Vratio is the ratio of the wind velocity at 
infinity to the escape velocity, expected to be in the range 
1-3. I is the integral along the line of sight to the pulsar 
from infinity, given by 

1=1 ^ 4^^. (2) 

with 7?2 the radius of the companion and r the distance 
from the centre of mass of the star. 

^ See |http://tempo.sourceforge.net] 
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We use this model to investigate offsets in DM over 
the entire orbit of the binary system given orbital param- 
eters determined through timing and integrating numeri- 
cally through the wind along the line of sight to the pul- 
sar. 

We calculate DM offsets by looking at four con- 
secutive orbits (MJD 52046 to 52971) that have data 
spanning multiple well-separated frequencies and bin- 
ning the TOAs into groups with multiple frequencies 
where possible, with lengths of five days or fewer. We 
run TEMPO with DM arbitrarily fixed at 738.5 cm~^ pc, 
allowing DM offsets to vary independently in each 
bin. The TEMPO error bars are doubled, and, using the 
PyMC Markov C hain Monte Carlo (MCMC) algorithm 
jPatil et alj|2010l) . we fit these offsets to the wind-model 
predictions. The orbital inclination angle is varied as a 
free parameter in order to find a best fit. From the incli- 
nation angle and the assumption of a 1 .4 Mq pulsar, the 
companion's mass is calculated using the binary mass 
function, and this in turn is used to estimate the compan- 
ion's radius, i?2. This is done by interpolating over the 
relevant ran ge of main-sequence masses in Table 15.8 of 
ICwJ ( I2OOOI) . Along with the inclination angle, we vary 
the quantity M /vrauo in Equation[T](equivalent to scal- 
ing the height of the curve in Figure [TJ and the base- 
line of the DM offsets (shifting the vertical position of 
the curve). The results of this fit are shown in the first 
two panels of Figure [T] with best fits of 53 ± 7 degrees, 
Qlg X 10"^° Mq yr"\ and 0.4ll2;}| cm"^ pc for the 
orbital inclination, the quan tity M /vratio, and th e DM 
baseline offset, respectively. iDemorest et al] i2010l) have 
measured a neutron star mass of 2 Mq, and so we also 
ran the MCMC code assuming a pulsar mass of 2.0 Mq. 
All fit parameters were consistent within the error bars 
of the results listed above, and so we continue to assume 
a pulsar mass of 1.4 Mq. 

It is apparent from Figure [T] that the model and ob- 
servations agree well. The fit has a reduced of 1.04. 
Using all three output fit parameters, we numerically cal- 
culate a predicted DM offset for each TOA. These off- 
sets are applied to all TOAs in TEMPO to remove the ef- 
fects of a varying DM from the timing residuals, even for 
epochs with only single-frequency observations. 

We can use the fit parameters to estimate some 
further quantities. Using the inclination angle obtained 
through this fit, and assuming a pulsar mass of 1.4 Mq, 
we estimate the mass of the companion to be 20I4 ''^o- 
This is con sistent with t he ear ly main-sequence star hy- 
pothesis of IStairs et al.l 1 200 ih. and with t he isochrones 
plotted in Figure 3 of iBassa et al.l ( I2OIII) for a main- 
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Figure 1. Top two panels show predicted dispersion measure 
offsets fit to measured changes in DM over four orbits of the 
J1740— 3052 system, both as a function of MJD and orbital 
phase. The orbital inclination determined by this fit is 53 ± 7 
degrees. The offsets shown are in reference to an arbitrary 'zero' 
value of 738.5 cm"'' pc. The upper horizontal axis on the top 
plot denotes January 1 each year The bottom panel shows aver- 
age scattering timescales of BCPM data observed at 1190 MHz 
varying across orbital phase. The data have been binned so that 
each point represents five measurements, except for the final 
point before phase zero, which represents four Note the increase 
in scattering neai' the phase of maximum DM, at which we ex- 
pect the beam of the pulsar to traverse more of the companion's 
wind than elsewhere along the orbit. In all panels, dotted vertical 
lines denote periastron. 



sequence companion younger than 10 Myr, which fits 
with the pulsar's characteristic age of 0.35 Myr 

Allowing a VraUo value in the range 1-3, our mea- 
surement of M / Vratio gives a mass-loss rate in the range 
(0.6-3)xl0~® Mq yr~^ for the companion, assuming a 
completely ionized wind. While this is lower than mod- 
els predict for radiative ly-driven winds by at least an 
order of magnitude (e.g. IVink et al]|200(]| : iMuiires et alj 
I2OI2I) , some Galactic late O-type main-sequence stars 
have in recent years been observed to have unexpectedly 
low mass-loss rates in j ust th is range (e.g. iMotch et al] 
I2OO5I : iMarcolino etal] |2009|) . This is known as the 
'weak- wind problem'. The cause of these weak winds is 
not established, but models that initiate winds using line 
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acceleration in addition to radiative acceleration do ap- 
pear t o reproduce the weak -wind problem for late 0-type 
stars iMuiires etai]|2012h . Whatever the cause, unless 
the wind is only partially ionized, it appears as though 
the companion to PSR J1740— 3052 may experience this 
phenomenon. 

The baseline offset is best interpreted in conjunction 
with the DM measured by TEMPO, and so we withhold 
discussion of this parameter until the next section. 

Finally, we find evidence for orbital variations not 
only in DM, but in the scattering timescale of pulses. In 
the third panel of Figure [T] we plot scattering timescales 
for 11 90 MHz data as a function of orbital phase, binned 
to increase signal-to-noise. It is clear that the amount of 
scatter increases where the pulses must traverse more of 
the stellar wind, which is the same phase at which the 
DM variations peak in the previous panel. 



2.5 Timing results 

We run T EMPO using the pulsar-main-sequence binary 
model of lWexl h998l) . The parameters fixed and fit in 
TEMPO are shown in Table [2] and the residuals for ob- 
serving frequencies greater than 1 GHz are plotted in 
Figure |2] While the fit includes ten spin period deriva- 
tives (with higher-order derivatives fit as a means to re- 
duce low-frequency timing noise), only the first two are 
shown. 

The lower-frequency data do not appear to follow 
the variation in arrival times between data at differ- 
ent frequencies as the > 1 GHz data do, and so our tim- 
ing solution excludes these data. Residuals that include 
these points are shown in Figure [3] to illustrate the poor 
fit. These TOAs account for timing offsets due to their 
large scattering timescales, as described in section lZ2l 

FITSCATTER is producing TOAs in these low- 
frequency bands that are late compared to the rest of 
our timing solution by up to ~7% of the pulse period. 
Profiles in this frequency range are very noisy, making 
it difficult to investigate the source of this delay. Per- 
haps the model used for scattering in the interstellar 
medium cannot be applied directly to scattering caused 
by a stellar wind. Some pulsars also show significant 



variations in profile shape between widely-spaced fre- 
quenc ies (e.g. iHassall et alj|2012l : IPhillips & Wolszczad 
this could affect timing at lower frequencies if 
present and undetected. 

The DM value output by TEMPO represents the 
baseline from which our orbital-phase-dependent DM 
offsets are measured, and this allows us to estimate how 



Table 2. Measured and derived TEMPO fit parameters, with la 
en'ors reported by TEMPO doubled. Ten spin period derivatives 
were used in the fit, but only the first two are reported here. 
Numbers in parentheses refer to uncertainty in the last reported 
digits. 

Measured Parameters 



1) 



RA (J2000)" 
Dec. (J2000)'" 
Period, P (s) 
yet Period Derivative, P 
2nd Period Derivative, P (s~ 
Epoch of Period (MJD) 
Dispersion Measure (cm~^ pc)** 
Binary Model 
Orbital Period, Pf, (d) 
Projected Semimajor Axis, x (Is) 
Eccentricity, e 

Longitude of Periastron, uj (° ) 
Epoch of Periastron, To (MJD) 
Advance of Periastron, ui (°/yr) 
Derivative of Projected 

Semimajor Axis, x 
Second Derivative of Projected 

Semimajor Axis, x (s"^)"^ 
Derivative of Eccentricity, e (s~^)'^ 
Derivative of Orbital Period, Pf, 
Second Derivative of 

Periastron, ui (rad s~^)° 
Data Span (MJD) 



17'"40™50!001 
-30°52'04'.'3 
0.570313411724(3) 
2.5504275(95) x 10"" 
9.1(6) X 10^26 
53191.0 
738.73(8) 
Main Sequence Star 
231.029630(2) 
756.90794(14) 
0.57887011(19) 
178.646811(17) 
52970.719801(11) 
0.000112(6) 

9(1) X 10-12 

< 3 X 10^20 

< 4 X 10-15 

< 3 X 10-3 

< 4 X 10-23 
50760 - 55622 



Derived Parameters 



Mass Function (Mq) 

Inclination Angle, i (°) 

Mass of companion, m2 (Mq ), 

assuming mi = 1.4 Mq 
ISM-only Dispersion Measure 

(cm^3 pc) 



;.722676(3) 
53(7) 



20 



+6 

4 



739.1(2) 



" Position held fixed in TEMPO; values from lBassa et al.l i201 ih 
^ The DM value output by TEMPO is somewhat arbitrary 
in our case, as discussed in the text. 

Fit while holding all other parameters fixed at the values 
shown. 

much of the DM is due to the interstellar medium (ISM) 
rather than the companion's wind. We expect it to be 
close to the arbitrary value of 738.5 cm-"^ pc we used 
as a baseline when determining the offsets, but there is 
no reason it should be the same, as we are now apply- 
ing a model DM offset to our entire set of TOAs. Indeed, 
we measure 738.73(8) cm""^ pc. As stated in the previ- 
ous section, the theoretical curve of Figure [T] was shifted 



up by 0.41 



+ 0.15 



cm pc to fit the measured DM off- 
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1300-1550 MHz 
1550-2000 MHz 
> 2000 MHz 




55000 



Figure 2. Best-fit residuals using DM-corrected data observed at frequencies greater than 1 GHz. The upper horizontal axis denotes 
January 1 each year, and vertical dashed lines show times at which periastron occurs. Note the quasi-periodic signal that repeats over 
slightly more than two orbits of the binary system. This is discussed further in the text. 



sets, and so our offsets are too high by this amount, or 
equivalently, our baseline DM value is too low by this 
amount. Adding this value, therefore, to the DM out- 
put by TEMPO, we obtain a value of 739.1(2) cm~^ pc, 
which is an estimate of the DM resulting only from 
the ISM. It should be noted that the companion wind's 
contribution should be included when observing PSR 
J1740— 3052, ideally with the orbital phase at the time 
of observation taken into account to dedisperse as accu- 
rately as possible. 

It is clear from Figure [2] that there remain some 
unmodelled systematics in our timing residuals, partic- 
ularly a prominent quasi-periodic signal that is present 
in data acquired at various frequencies and observato- 
ries. Figure |4] shows a segment of the residuals for a fit 
using only four spin derivatives, rather than ten, to em- 
phasize that we do not think these remaining systemat- 
ics are merely an artefact caused by fitting out a high- 
order polynomial. Using a Lom b-Scargle periodogram 
dScargldl 19821: IPress et alfcool § 13.8), we measure the 
period of this signal to be I.40lo;io years, or 2.22to\\ 
orbital periods, where the error bars are taken as the 
change in frequency (subsequently converted to period) 
required to reach half the peak value on the periodogram. 
This signal will be discussed further below. 



3 DISCUSSION 

3.1 Secular changes in orbital parameters 

The value measured for the advance of periastron 
here, 6j = 0?000112(6) yr-\ is si naller than 
the va lue 0?00021(7) yr"^ measured by IStairs et alj 
( l200lh . though they are consistent within twice the lat- 
ter's la error bars. Following the calculation in Equa- 
tion 2 of that paper, but with mi = 1.4 M0 and 
m2 = 20 M0, general relativity predicts an advance 
of [j = 0? 00026 yr~^. This is approximately double 
the value measured here, suggesting that there is a coun- 
teracting classical contribution due to a mass quadrupole. 
Conversely, if we assume our measured io to be due en- 
tirely to general relativity, we can determine the total 
mass of the system using this value along with the mea- 
sured orbital period and eccentricity. This provides a 3(t 
upper limit of 7.3 Mq, which is considerably smaller 
than the 8.7 M© binary mass function. Thus, without a 
counteracting classical contribution, we would require a 
neg ative pulsar mass. 

IStairs et al. I l l200lh show that for a B-type main 
sequence star companion, a mass quadrupole due 
to tidal effects is negligible compared to a spin- 
induced quadrupole. Following their methods, we find 
that the classical contribution to uj is approximately 
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Figure 3. Best-fit residuals using DM-corrected data observed 
at frequencies less than 1 GHz, shown as triangles, with verti- 
cal dashed lines showing times at which periastron occurs. The 
other points plotted are a segment of those in Figure |2] also 
corrected for DM variations. The low-frequency TOAs are pro- 
duced by FITSCATTER to account for scattering due to stellar 
wind and the interstellar medium. The amount of delay in these 
low-frequency residuals is considerably more than we would ex- 
pect even from our DM-variation model. 



—0? 00015 yr^^. Using this and our measurement of the 
rate of change of the projected semi-major axis x, for 
which there was previously only an upper limit, we find 
a greater angle between spin and orbital angular momen- 
tum vectors in the companion than was previously es- 
ti mated . As s hown in Figure |5] which uses the model 
of lWexl ( Il998h . the companion's spin momentum differs 
by 40°-87° from being either aligned or anti-aligned 
with the orbital angular momentum of the binary sys- 
tem. Such a large relative tilt would likely have been in- 
duced by an asymmetry or 'k ick' in the superno va ex- 
plosion that created the pulsar ( iKasDi etal.ll996iJ) . Mak- 
ing the assumption that the angular momenta in the pre- 
supe mova system wer e aligned and that the orbit was cir- 
cular jKaIogera|[T99^ . we find the kick velocity was at 
least 50 km/s. 

A known alternative mechanism for producing sec- 
ular changes in x and ijj arises simply from the change in 
view ing geometry du e to the proper motion of the system 
(e.g. lKopeikinlll996h . Proper motion has not been mea- 
sured in our system, but we can easily see that this effect 
may be disregarded. A transverse speed of 500 km/s, ab- 
surdly high for such a massive system, leads to a motion 
of about 10 mas/yr at the estimated distance of 11 kpc. 
This would contribute at most ~ 4 x 10~^ deg yr~^ to 
di and 9 X 10^^'^ to x. Comparing with our measured 
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\ Green Bank 
I Parkes 
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Figure 4. A segment of the residuals after fitting only four fre- 
quency deiivatives rather than ten, with colours representing dif- 
ferent telescopes. It is clear that the 2.2Pi, signal is present at all 
sites, and is not merely an artefact caused by fitting out a higher- 
order polynomial. Axis labels are the same as in Figure |2] and 
vertical dashed lines show times at which periastron occurs. 



values in Table |2] or with the classical contribution to w 
above, we see that this effect falls short of our measure- 
ments by an order of magnitude or more in each case. 
Thus, we consider it no further. 

The distribution of angles shown in Figure |5] leads 
to a 3a lower limit of 3.2 x 10~^ AU'^ on the quadrupole 
moment q of the companion. This is determined by al- 
lowing (jj, X, and the inclination angle to vary by 3cr and 
finding the allowed combination of angles 6 and <3>o that 
minimizes q in the equations of Wex ( 199 8). This cannot 
place an upper limit on q, but a Monte Carlo sampling 
of the parameter space gives a median value that is only 
about twice the lower limit at 6.7 x lO"'^ AU^, assum- 
ing a uniform distribution over the orbital-plane preces- 
sion angle <I>o. Taking this median g, companion mass 
20 M0 and a stellar radius of 8 Rq determined by inter- 
polati ng a stand ard table of main-sequence mass es and 
radii (|Coxll2000h . we can use Equation 6 of Stairs et al] 
( l200l[) to estimate roughly the stellar spin period Ps . We 
find Ps « 2 X 1 0^ s, which remain s consistent with the 
value assumed bv lStairsetai] ( l200lh . 



3.2 Low-frequency residual signal 

The source of the 2.2Pi, signal in the timing residu- 
als is not known. For the sake of interest, we imag- 
ine a third body orbiting the pulsar-main-sequence bi- 
nary system with this period. As a simple approxima- 
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Figure 5. Constraints on the angle of precession of the orbital 
plane, and the angle between the main-sequence star's spin 
and orbital angular momenta, d, given the masses of the two 
stars and the inclination angle of the system. The grey regions 
are excluded because x is positive and the hatched regions are 
excluded because the classical contribution to u) is negative. Be- 
cause the inclination angle is by far the dominant source of error, 
the solid lines are drawn assuming the nominal measured values 
of X and u). If the spin and orbital angular momenta are closer to 
being aligned than anti-aligned, we find that 40° < d < 87°, 
with 9 close to 80° for most values of <I>o- Otherwise, we may 
subtract these angles from 180° . After Wex 1998, Figure 1 1 . 



tion, we treat the binary system as a single body of mass 
20 + 1.4 = 21.4 Mq located at its centre of mass in 
a circular orbit with a planet. We also assume the incli- 
nation of this orbit to be the same as that for the binary 
system, which we fixed at 53°. We find that in such an 
approximation, given the uncertainties in inclination an- 
gle, main-sequence star mass, and the amplitude of the 
signal (which we estimate to be anywhere from slightly 
less than 1 ms to 2 ms), a planet would be 2-8 Earth 
masses, orbiting at a radius of 3^ AU. 

Given the inclination angle, the semimajor axis of 
the inner binary is 1.9 ± 0.2 AU, with a maximum sepa- 
ration of 3.0±0.3 AU. A circumbinary planet thus seems 
unlikely at such a short di stance, making this rough 
approximation unappealing. IPoolin & BlundelH ( 1201 ih 
have run simulations to test the stability of circumbinary 
orbits. The results of those simulations suggest that a sys- 
tem with the eccentricity and mass ratio of J1740— 3052 
could not host a stable circumbinary orbit closer than 
roughly 5 times the semimajor axis of the binary (e.g. 
Figure 14 of that paper). 

Furthermore, a planet orbiting only the pulsar with 



this period would be at roughly the same distance as the 
companion star, which is untenable. 

A more satisfying explanation would be one that 
does not invoke an additional orbiting body. If it were 
some consequence of the companion's spin, we would 
expect a much shorter period, as main-sequence sta rs ro- 
tate w ith periods on the order of days, not years (e.g. lCoxl 
l2000l. Table 15.8). In particular, in Section [3T| above. we 
estimated the spin period of the companion to be only a 
few days. 

lApplegate & Shahaiiil h994h explain such a 'super- 
orbital' period in eclipsing pulsar binary PSR B 1957+20 
by variations in the quadrupole moment induced by tidal 
activity. The far greater separation and companion mass 
in the PSR J1740— 3052 system make this seem a very 
unlikely mechanism for the observed variations, as well 
as our assumption in the previous section that tidal ef- 
fects are small in this system. The orbit of the pulsar be- 
ing tilted at least 40° and perhaps close to 90° from the 
spin of the companion would further reduce the effec- 
tiveness of such a mechanism. 

If the comp anion is a Be star, as in the case of 
PSR B 1259-63 ( Johnston et al.' 1992), its circumstellar 
disk may be tilted due to the misalignment of the star's 
spin and orbital angul ar momenta, and mi ght thus be ex- 
pected to precess (e.g. iMartin et al.ll201 ih . A precessmg 
disk would introduce a new source of periodicity into 
the system that mi ght explain the observed variations. 
iMartin et al] ( 1201 ih suggest that the B-star binary PSR 
J0045— 7319 should show oscillations on a timescale of 
about a year if a disk is periodically ejected, and the same 
might apply to J1740— 3052. Evidence of a disk around 
the companion star (such as hydrogen emission lines) 
would put constraints on such a scenario in our system. 



4 CONCLUSION 

We have updated the IStairs et al] ( 1200 ih timing model 
for PSR J1740— 3052 using data acquired on several in- 
struments. From the timing results, we see orbital-phase 
variations in the scattering and dispersion of the pulsar's 
radio beam that are clear indicators of a stellar wind com- 
ing from its massive companion, lending further support 
to the pulsar-main-sequence binary picture for this sys- 
tem. Furthermore, the dispersion-measure variations fit 
very well to a simple radiative stellar wind model for the 
charged-particle column density. The mass-loss rate de- 
rived from this fit is lower than predicted by models, but 
matches the rates observed in some late O-type main- 
sequence stars found in the Galaxy in recent years. This 
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weak-wind problem remains an area of active current re- 
search. If the companion of PSR J1740— 3052 is indeed 
such a star, then its spectral class is presumably late O- 
type rather than early B-type. 

With our new timing solution, we are able to mea- 
sure the derivative of the projected semi-major axis x, 
which leads to new constraints on the geometry of the 
binary system — the pulsar's orbit appears to be tilted 
substantially from being aligned (or anti-aligned) with 
the companion's spin. This tilt of the orbital plane was 
likely caused by asymmetries in the supernova explo- 
sion that produced the pulsar, implying a kick of at least 
~50km/s. 

The origin of the quasi-periodic 2.2Pt si gnal seen in 
the tim ing residuals of Figure[2]is a mystery. Bassa et ahl 
( I2OI ih suggest that spectroscopic observations of the 
companion could be made with the aid of adaptive op- 
tics, further constraining the orbital parameters and spec- 
tral class of the star. Such constraints might narrow down 
the possible mechanisms for producing this signal, par- 
ticularly if evidence is found for a circumstellar disk. 

Finally, considering the winds of both the pul- 
sar and its main-sequence companion, we might ex- 
pect a shock front between the two to produce emis- 
sion at X-ray and gamma-ray energies in this bi- 
nary system. High-energy emission has been observed 
in associati on with the pul sar-main-sequence binary 
B 1259— 63 dTavani et al . 1994) and the millisecond pul- 
sar B 1957+20 jStappers et al. .2003). in both case s at- 
tributed to a shocked pulsar wind. Stairs et al] ( I2OOII) ex- 
amined archival X-ray observations in the field of PSR 
J1740— 3052 and found no significant emission, but were 
unable to rule out the existence of shock emission. 

PSR J1740— 3052 is one of just a few systems 
that represents a particular early stage in the evolution 
of binary systems containing neutron stars. Most such 
binaries presumably begin as a pair of massive main- 
sequence stars and pass through a relatively brief phase 
during which one has evolved to form a neutron star, 
while the other still burns hydrogen in its core. Not only 
does this strengthen our models of binary evolution; as 
seen here, it provides a unique means of probing the wind 
and gravitational properties of a main-sequence star. 
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